Statistical mechanics of double-stranded semi-flexible polymers 
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We study the statistical mechanics of double-stranded semi-flexible polymers using both ana- 
lytical techniques and simulation. We find a transition at some finite temperature, from a type of 
short range order to a fundamentally different sort of short range order. In the high temperature 
regime, the 2-point correlation functions of the object are identical to worm-like chains, while at low 
temperatures they are different due to a twist structure. In the low temperature phase, the polymers 
develop a kink-rod structure which could clarify some recent puzzling experiments on actin. 
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A double-stranded semi-flexible polymer chain is the 
basic structure of many biopolymers of which examples 
are DNA and proteins such as actin. There has been a lot 
of recent experimental interest in systems of biopolymers 
ranging from the elasticty of biopolymer networks and 
its use in the prediction of the mechanical properties of 
cells to the direct visualisation of single chain properties. 
The model most used in the study of biopolymers is that 
of the worm-like chain jjj in which the polymer flexibil- 
ity (structure) is determined by a single length, the per- 
sistence length £ p which measures the tangent-tangent 
correlations. For example, DNA has a persistence length 
ip ~ 50nm whilst for actin £ p rs 17/im. These biopoly- 
mers are known to have a more complex 'twisted' struc- 
ture. The multi-stranded nature of these polymers is also 
not taken into account in a simple worm-like chain model, 
ft is not clear if such a fine structure will have an effect 
on the global properties of these objects. A possible ef- 
fect of such fine structure is what we attempt to study 
in this Letter. Our model is, in a sense, microscopic be- 
cause the interaction between the bend and twist degrees 
of freedom is a result. This is fundamentally diffent from 
previous approaches 0] which try to include the twist 
degrees of freedom by adding extra terms to the free en- 
ergy. 

We study a version of the rail-way track model of 
Everaers-Bundschuh-Kremer (EBK) j|] for a double- 
stranded semi-flexible polymer, embedded in a d- 
dimensional space for arbitrary d. Excluded volume and 
electrostatic interactions have been ignored. We find that 
the system, has qualitatively different properties in the 
low temperature and high temperature regimes, in con- 
trast to what one might naively expect from an inherently 
one-dimensional system with local interactions and con- 
straints. The tangent-tangent correlation function de- 
cays exponentially in the whole range of temperatures 
with a "tangent-persistence length" £tp that has a very 
slow temperature dependence, and whose scale is deter- 
mined by the (bare) persistence length of a single strand 
£ p = n/ksT (k is the bending stiffness of a single strand). 



Note that it is independent of a, the separation of the 
two strands, which is the other relevant length scale in 
the problem. However, the correlation function of the 
bond-director field, defined as a vector that determines 
the separation and coupling of the two strands of the 
combined polymer system, has different behaviour below 
and above the temperature T c ~ 4.27n/dk B a. While it 
decays purely exponentially for T > T C1 there are addi- 
tional oscillatory modulations for T < T c . The related 
"bond-persistence length" £qp does not change appre- 
ciably at high temperatures , where its scale is again 
set by £ p alone. In the low temperature phase, how- 
ever, £bp does show a temperature dependence. In par- 
ticular, £ BP ~ 4 /3 ° 2/3 <x 7 1 - 1 / 3 for T -> 0, while 
£bp ~ £ P for T ~ T c . Similarly, the "pitch" H, de- 
fined as the period of oscillations in the low temperature 
regime, changes drastically with temperature , ranging 
from H ~ £ B p ~ £l /3 a 2 / 3 near T = 0, to H ~ near 
T = T c . At T — we regain a flat ribbon which has true 
long-range order in both the tangent and bond-director 
fields. The ribbon is essentially a rigid rod. As we ap- 
proach T = 0, the persistence lengths and the pitch di- 

~l/3 

verge with the scaling H ~ £bp ~ t T ' p . 

The system is composed of two semi-flexible chains, 
each with rigidity k, whose embeddings in d-dimensional 
space are defined by Vi(s) and r 2 (s). The Hamiltonian 
of the system can be written as the sum of the Hamil- 
tonians of two worm-like chains subject to the inexten- 
siblity constraints [0. The ribbon structure is enforced 
by having ra(s') separated from ri(s) by a distance a, 
i.e. T2(s') = ri (s) + an(s) where \s— s'\ can be non-zero 
but is small. We have defined a bond-director field n(s), 
a unit vector perpendicular to both strands. The chains 
are assumed to have "permanent" bonds (e.g. hydrogen 
bonds) that are strong enough to keep the distance be- 
tween the two strands constant. In Ref. H, it is argued 
that the relevant constraint on the system would then 
be that in a bent configuration the arc-length mismatch 
between the two strands should be very small. We can 
calculate the arc-length mismatch for the bent configu- 
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ration as As = |r2(s) — ri(s) + a n(s)|, where a is the 
separation of the strands. We impose the constraint as a 
hard one, namely, we set As = 0, as opposed to Ref. j|. 
Physically this means we do not allow bends in the plane 
of the ribbon. These bends are less important in d > 2 
because as we shall see the lower length-scale will be set 
by the 'pitch' which will make the in plane fluctuations 
of the ribbon irrelevant This simplifying assump- 
tion should not change the behavior of the system ||. 
We implement the constraint As = by introducing the 
"midcurve" r(s): ri(s) =r(s) + | n, r 2 (s) =r(s) — | n. 
In terms of the tangent to the midcurve t = dr/ds and 
the bond-director n the Hamiltonian of the system can 
now be written as 



n = 



dt(s] 
ds 



a 

T 



d 2 n(s) 
ds 2 



(1) 



subject to the exact (local) constraints 



a dn, 9 



1. 



1, (t±^)-n = 0. (2) 



This completes the formulation of the model. 

The statistical mechanics of semi-flexible chains are 
difficult due to the constraint of in-extensibility. Vari- 
ous approximation methods have been devised to tackle 
the problem. A successful scheme, that somehow man- 
ages to capture the crucial features of the problem, is 
to impose global (average) constraints rather than lo- 
cal (exact) ones g. This approximation is known to be 
good for calculating the average end-to-end length but 
not so good for the whole distribution. It corresponds 
to a saddle-point evaluation of the integrals over the La- 
grange multipliers, that are introduced to implement the 
constraints [Q. In this sense, it is known to be a "mean- 
field" approximation in spirit. To study the effects of 
fluctuations on the mean-field result, we have performed 
a l/d-expansion similar to the one successfully used by 
David and Guitter to study the crumpling transition of 
crystalline membranes ||. We see that no divergent be- 
haviour appears in the diagrams of the 2-point correla- 
tion functions which means that the mean-field behavior 
of these functions is not changed by fluctuations. This 
does not preclude differences in higher order correlation 
functions. 

With the above discussion as justification, we ap- 
ply the same approximation scheme to our prob- 
lem defined above. The local constraints in Eq.(||) 
are relaxed to global ones. To do this, we add 
the corresponding "mass terms" to our Hamiltonian 
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where 6, c, and e 



are dimensionless constants. We then determine the con- 
stants self consistently by demanding the constraints of 
Eq.(|]) to hold on average, where the thermal average is 
calculated by using the total Hamiltonian H + H m - Note 
that in choosing the above form, we have implemented 
the "label symmetry" of the chains, namely, that there 
is no difference between two chains. The self-consistency 
lead to the following set of equations for the constants b 



and 



a 2 V^ _ 



c(b + Vc) = 



and e = 0. 



4V2b ' idP p d ' "V 1 V -/ 2a 4 

The above equations, which are nonlinear and difficult 
to solve exactly, determine the behavior of b and c as 
a function of u = a/£ p . We have solved them numeri- 
cally in d = 3 and the solutions are given in Fig. 1. One 
can solve the equations analytically in two limiting cases. 
For u < 1 we find b = d 2 /32 , and c = (d/ V2) i/3 u~ s/3 , 
whereas for m> lwe find b — d 2 /8, and c = 4/u 4 . In 
Fig. 1, the behavior of b and c is plotted as a function of 
u. Note that u is proportional to T and can be viewed 
as a measure of temperature. 




u = all 



FIG. 1. The solution of the self consistent equations for 
the constants b and c as a function of u — a/lp in d = 3. 
The value u c ~ 4.27/d corresponds to the transition point. 

We can then calculate the correlation functions. For 
the tangent-tangent correlation one obtains 



(t(s)-t(0)) = 



a / j— t 

— = exp -V2& - 
4yM V 1 



(3) 



whereas for the bond-director field one obtains 



(n(s) • n(0)) 



dtl 



exp 



(-{b- ^/W~c) 1 / 2 exp (-(6+ VW^c) 1 / 2 



2a 2 Vb 2 



{b - V6 2_ ^) 1 /2 



(b + vv^zyi 2 



(4) 
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The tangent-tangent correlation (Eq.(||)) is exactly 
what we obtain for a single worm-like chain, and im- 
plies uniform behavior for all temperatures. Eq.(|^) on 
the other hand, indicates a change of behavior at b 2 = c 
for the bond-director correlation. The correlation is 
overdamped for b 2 > c (high temperatures), while it 
is underdamped (oscillatory) for b 2 < c (low tempera- 
tures). The interesting point b 2 — c happens for u c — 
2 19 / 4 /(2 + V2) 3/2 d ~ 4.27/rf, that leads to the value for 
T c quoted above (see Fig. 1). We also find a divergence 
in the specific heat, Cy — where F — — ksT log Z 
at T c . It should be noted that it is not a thermodynamic 
phase transition in the sense of long-range ordering and 
broken symmetry. It is a cross-over that appears due to 
competing effects, and the transition is from a state with 
some short-range order to a state with a different short- 
range order. Similar phenomena have been observed in 
Ising-like spin system with competing interactions j?j and 
the cross-over (transition) point corresponds to a type of 
'Lifshitz point ' for a 1-d system. 

The nature of competition in our double-stranded 
polymer system can be understood using a plaquette 
model. Consider rectangular plaquettes attached to each 
other to form a ribbon. The attached sides (edges) cor- 
respond to the bond-director field, while the other sides 
determine the tangent-director field. The twist degree of 
freedom corresponds to twisting the plaquettes against 
each other. The energy expression corresponding to 
bends comes from the product of the tangent-directors 
of the neighbouring plaquettes, and has no competition. 
On the other hand, the energy for twists that comes from 
the product of the bond-director fields of the neighbor- 
ing plaquettes, does show competition. These competing 
effects come from an effective interaction between bond- 
directors that are next nearest neighbors. This is due 
to the fact that alike twists meeting at an edge tend to 
unwind (annihilate) each other, while unlike twists when 
they meet arc trapped; they do not annihilate each other. 
This competition is only present at nonzero temperatures 
and is merely due to topological constraints of the ribbon. 

It is useful to study the bond-director correla- 
tion in the limiting case b 2 <C c, that corresponds 
to relatively low temperatures. Using the asymp- 
totic forms for b and c, one obtains (n(s) • n(0)) = 

^2 exp^( 3 ^) 1/3 S ) sin^^^ + f), for 

very low temperatures. From the above expressions for 
the correlation functions, one can read off the persistence 
lengths £tp and £bp, and the pitch H, as summarized 
above. 

An intriguing feature of the behaviour of this model is 
that, although the ground state (T — 0) configuration of 
the system is a flat ribbon, and supports no twists, upon 
raising the temperature, a twisted structure with short 
range order emerges. We have confirmed this by perform- 



ing extensive Molecular Dynamics (MD) /Monte Carlo 
(MC) simulations of double-stranded semi-flexible poly- 
mers. A bead-spring model with bending and stretching 
energies was used. We combined a velocity- Verlet MD 
coupled to a heat bath with an off-lattice pivot MC algo- 
rithm. The MD was useful for equilibriating the shorter 
length-scales and MC the long length-scales. Details of 
the algorithm and computational method will be pre- 
sented elsewhere 0. 




FIG. 2. Typical conformations from MD/MC simulations 
of a ribbon made up of two chains of 400 monomers (a) above 
(b) near and (c) below T c 

We show typical cquilibriated conformations in Fig. 2. 
The snapshots of the polymer configurations suggest that 
at low temperatures the polymer can be viewed as a col- 
lection of long, twisted (straight) rods that are connected 
by short, highly curved sections of chain which we call 
"kinks" . This picture can be accounted for using a simple 
argument. We can model our system of two semiflexible 
polymers subject to the constraint of constant separation, 
as a semiflexible ribbon, i.e. a semiflexible linear object 
with anisotropic rigidities whose Hamiltonian reads 




where Kij = «|| ni,nj+Kj_ (dij—titj—niTij) determines the 
rigidity anisotropy of the ribbon, corresponding to bend- 
ing parallel or perpendicular to the bond-director field. 
The ribbon structure would require km k±. To be con- 
sistent with the hard constraint (see above) of constant 
separation of the polymers, we should take the limit of 
infinite /cm. The Boltzman weight with the above Hamil- 
tonian in the km — > oo limit can be easily shown to require 
the constraint d v a ^ • n(s) = 0, to hold exactly at every 
point of the ribbon. Recalling dt/ds = H(s) e from the 
Frenet-Seret equations, where H(s) is the curvature at 
each point and e is a unit normal vector to the curve, 
we can write the constraint as H(s) e(s) • n(s) = 0. 
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It tells us that at each point either H (s) = 0, that is 
we have straight, twisted sections with no curvature, 
or e(s) • n(s) = 0, curved regions where the the bond- 
director is perpendicular to the curve normal. The for- 
mer case would correspond to the rod segments whereas 
the latter would correspond to the "kinks" . Performing 
the configuration sum we obtain that the average sep- 
aration of the kinks is equal to the tangent-persistence 
length, which is quite large at small enough tempera- 
tures. As the temperature is raised, the kinks get closer 
so that at some temperature the average kink separation 
becomes comparable to the size of the kinks, where the 
rod-kink pattern disappears. This analysis can be under- 
stood in the context of the mean-field (n, t) model above 
by observing that at low temperatures £bp <C £tp, one 
can imagine that there are roughly speaking rod-like seg- 
ments of length £tp , each supporting a number of shorter 
segments of length ^bp that are twisted, but de-correlated 
with one other. As the temperature is raised, the num- 
ber of twisted rods in each segment N — ^tp/^bp de- 
creases very quickly until it saturates to unity at T = T c . 
For higher temperatures the mechanism changes, and the 
bond correlations are cut off by the tangent fluctuations. 
Hence, the short-range twist order does not exist any- 
more. All the main features of the above picture have 
been observed in the simulation. We plot the (n(s)-n(O)) 
correlation function from the simulation in Fig. 3. 



sive MD/MC simulations. We have shown novel non- 
trivial differences between the high, low and zero temper- 
ature behaviour. An alternative formulation of the prob- 
lem (see Eq.(||) and the following paragraph) where the 
constraints are implemented in a different way also sug- 
gests qualitatively similar results [||. We also note that 
the low temperature phase can be used to explain some 
recent experiments done on single actin filaments by Kas 
et al I! who observed a length-scale dependent bending 
rigidity as well as the kink-rod structure explained above. 
The high-q small bending rigidity can be associated with 
the kinks and the low-g large bending rigidity with the 
rods. Of course, upon averaging the object would have 
worm-like chain statistics. Finally we mention that there 
are lots of interesting prospects for the study of such com- 
plex linear objects. It is known that biopolymers, such as 
microtubules, can also be multi-stranded objects. This 
method could be easily be extended to describe triple- 
stranded objects. The effect of an intrinsic twist changes 
the ground state but does not change any of the conclu- 
sions of our description |8| though we expect it to make 
the effective persistence length much higher. 

We have benefited from numerous discussions with R. 
Bundschuh, R. Ejtehadi, R Evaerers, E. Frey, G. Grest, 
M. Kardar and M. Piitz. We would like to thank the 
hospitality of ICTP, Trieste, where some of this work 
was done. 
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FIG. 3. The (n(s)-n(O)} correlation function measured 
in the simulations for temperatures /3 = 1/ksT — 1,5, 50 
corresponding to b 2 > c, b 2 « c and b 2 < c. The averages 
were done over ~ 10 4 statistically independent samples. 

In conclusion, we have studied the properties of a well- 
defined model of a double-stranded semi-flexible polymer 
using a mean-field analytical approach as well as exten- 
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